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Abstract 



A detailed study of an S — % spin ladder model is given. The ladder 
consists of plaquettes formed by nearest neighbor rungs with all possible 
S{7(2)-invariant interactions. For properly chosen coupling constants, 
the model is shown to be integrable in the sense that the quantum Yang- 
Baxter equation holds and one has an infinite number of conserved 
quantities. The i?-matrix and L-operator associated with the model 
Hamiltonian are given in a limiting case. It is shown that after a simple 
transformation, the model can be solved via a Bethe ansatz. The phase 
diagram of the ground state is exactly derived using the Bethe ansatz 
equation. 
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Heisenberg spin ladders and generalized spin ladders have attracted considerable atten- 
tion in recent years, due to the developing experimental results on ladder materials and the 
hope to get some insight into the physics of metal-oxide superconductors [j]. Especially 
generalized ladders including other couplings beyond the simplest case of rung and leg ex- 
change interpolate among a variety of systems and exhibit a remarkably rich behavior [2-9] . 
In particular, it has been shown that the diagonal interactions may cause frustration and 
change the structure of the ground state |2|, [|, while the biquadratic interactions, which can 
arise due to effective spin-spin interaction mediated by phonons in real magnetic systems 
PJ, tend to produce dimerization and may lead to a phase transition into a "non-Haldane" 
spin liquid state with absence of magnon excitations J|, |[. 

As spin ladders are generally not equivalent to spin chains with nearest neighbor in- 
teractions, till now little is known about integrable spin ladder models^. In this letter we 
study a generalized S — ^ spin ladder system with both isotropic exchange interactions and 



biquadratic interactions. Using ideas related to the quantum Yang-Baxter equations [10 



we found in our systems some cases of integrable ladder systems, in the sense of models 
having an infinite number of conserved quantities with explicit R matrices satisfying the 
Yang-Baxter equation. Properly choosing the spectral parameter, we get a Hamiltonian 
consisting of only nearest-neighbor and next-nearest-neighbor interactions. This model can 
be solved via an ordinary Bethe ansatz. 



4 The present work was submitted for publication in Euro. Phys. Lett.. During the preparation of a 
revised version of this paper, another integrable ladder model without diagonal interactions was presented 
in fill. 
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We consider a symmetric 16 x 16 matrix: 
R(x) = 
( a\ 

a 3 9ez2 3a2 3^2 
9a2 03 3&2 3a2 

fli 

3a 2 36 2 05 02 

I62 4a 2 46 2 2a 2 

4a 2 a 4 8a 2 46 2 

a 5 a 2 3a 2 36 2 

3b 2 3a 2 a 2 a 5 

46 2 802 04 4a 2 

2a,2 46 2 4a2 16x 

a 2 a 5 36 2 3a 2 

«i 

3a 2 36 2 «3 9a 2 

36 2 3a 2 9a 2 a 3 

V a i / 

where ai = 2(— 1 + 9x), a 2 = —b 2 = (—1 + x), a 3 = 7 + 9x, a 4 = 2(3 + 5x), a 5 = — 1 + 17x 
and a: G C. 

Let denote a 4-dimensional complex vector space. R(x) takes values in End^iV ® V) 
and satisfies the quantum Yang-Baxter equation (T(| : 



R 12 (x)R 23 (xy)R 12 (y) = R 23 (y) R 12 (xy) R 23 (x) , (1) 

where Rij denotes the matrix on the complex vector space V ®V ®V , acting as R(x) on 
the i-th and the j-th components and as the identity on the other components. Namely 
R12 — R ® 1, -R23 = 1 ® R and 1 is the identity operator on V. 

Let us set Ti = Vi®V 2 ® ... ®Vn, N e IN. The corresponding L-operator acting on the 
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) 
) 

)21 = 
) 



i-th space Vi, % — 1, 2, N, is a 4 x 4 matrix with entries 

(A(x)) n = |(4(3x-l)a^ + ci(^^ + ^i>2°.) 

+66 2 (<7 1 - i <7+ + (7+ <7 2 - ) + 2(1 + 3x)a* lti aZj 
|(a 5 ((7^(7 2 - + cr^o-^) + a 2 (^4i + ^M* 7 !,*)) 
|(a 2 ((7^(7 2 - + crf^) + a 5 (cr 1 - <7^ + cr^ af^)) 
2ci(j^(j^ 

|(a 3 ((7^(7^ + al^d + 9a 2 (cr^(7^ + af^)) 
§((15* - 7)<7°^ + 4(1 + *)« i( 7 2 * - a*.**.) 
+8a 2 (cr 1 + - cr^ <7^) - CiO^ct^) 

i(36 2 (l + el^i) + 2a 4 a^a+ t + ^a 2 (af ti a^ + a\^ z %i - a\^) 

U 9a ^ a ii a t + °U a ti) + a ^tA,i + crjicrf,.)) 
5(362(1 + + ^(a^a^i + al^li - o\^ z 2 ^ + 2a 4 <7£ i <7^) 

±((15* - 7)^(7", + 4(1 + *)(<7°^ - af,^) 
+8a 2 (<7i; i <7j i - (7^(7^) - Ciaf ^l J 

|(a 3 ((7^(7 2 - - <7^<7 2 ~ ) + 902(0-^(7^ - (T^,)) 
2Cl(7^(7^ 

l(a 2 (a^a 2 + - o\^$^ + a 5 (<rj i <r§ )i - <rt^2,i)) 
±(a 5 K^ 2 + - (7^(7+ ) + a 2 {at^% - 
|(4(3a; - lKyo, - c^a^ + af^) 
+6fc 2 (<7 1 >+ + a+ <7 2 - ) + 2(1 + Zx)al^ 

where c\ = ai/2, af^ a z Qi (resp. £7$ J, 9 = 1,2, are Pauli matrices (resp. 2x2 identity 
matrix) acting on the space V$. A direct calculation shows that L^x) satisfies 
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(A 
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(A 
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J32 — 

)33 = 



34 
41 
42 
43 
44 



£ (J J [A(x) ® L t (y)} = My) ® A(*p (J J , 
where *, y G C, y 7^ 0, <E> is the tensor product of matrices. Let 



1,2, ...,7V, 



(2) 



T(x) = L N (x)L N ^ 1 (x)...L x (x). 



(3) 



We have the fundamental commutation relations given by: 

h(A [T{x)®T{y)] = [T{y)®T{x)]R^- 



(4) 



Let t(x) = TroT(x), where Ttq takes trace on the space of 4 x 4 matrix. According to 
(fD a system with Hamiltonian of the form H XQ = J-^ \ogt(x)\ x=xo for some x G € (such 
that \ogt{x) is defined and differentiable, e.g. x = 1, see @) and real constant J, acting 
on 7i, has an infinite number of conserved quantities t(x): 

[H xo ,t(x)]=0, VxeC, (5) 

where [ , ] stands for the commutator. A system with Hamiltonian H xo is then by definition 
an integrable system. For arbitrary value of Xq, H xq generally describes integrable models 
with long range interactions. These models can be exactly solved by using the algebraic 
Bethe Ansatz method. 

Let P denote the permutation matrix on the ladder and set R(x) = PR(x). We see that 
R(x)\ x= i is proportional to the permutation matrix P. Therefore for xq = 1, the Hamiltonian 
H = Hi describes a system with nearest-neighbor interactions: 

d N 

H = —logt(x)\ x=1 = y2Ki+i 
dx ^ 

JV 

= y^[5(Sx,j • S2,« + Si 5 i + i • S2,i+i) + 3(8^ • Si^+i + S2,j • S2,i+i) (6) 

i=l 

— 3(81^ • S2,j+1 + S2,j • Si^+i) — 12(S2,j • Si j j + i)(Si j j • S2,j+l) 

+20(S M ■ S 2li )(S lii+1 ■ S 2 , m ) + 12(S M ■ S M+1 )(S 2ji • S 2 , m ) + f 1 ® 1], 
where the periodic condition is assumed: Sg.jv+i = Se,i, = 1,2 and S^^ = (cr^, o V q i; <J z e J/2 
is the spin vector operator on V^. Taking Si^ (resp. S2,») to be the spin operator on the first 
(resp. second) leg of the z-th rung of a ladder, and to be the tensor space for the action 
of these two spin operators, the Hamiltonian (||) describes an integrable spin ladder system 
with periodic boundary conditions. The model is S'?7(2)-symmetric, i.e., [H, S 1 } = 0, where 

N 

S l = XX^M + ^2 i)> ^ = x i Vt z are ^ e total spin operators of the ladder. 



i=i 



Remark: Like the well known Aflieck-Kennedy-Lieb-Tasaki model ||12|| , the model (P) 
derived from the transfer matrix has no free parameters. Nevertheless, as we shall show 
below, it can be generalized to a model with two free parameters, i.e., the coupling constant 
along rungs and the coupling constant of the rung-rung biquadratic interactions without 
losing the integrability. The present model and the ones discussed in the previous papers 
||, |9| have same interaction terms in Hamiltonian but with different coupling constants. The 
models in [5,9] with more free parameters are generally not integrable. 



To exactly solve the model (|6]) the analytic algebraic Bethe Ansatz method may be ap- 
plied. The reference state with all the spins up is an eigenstate. Some degenerate eigenstates 
states can be obtained by applying the operator S~ = S x — iS y . The combinations of the 
products of (T(x))i2 and (T(x)) 13 can be used to construct "Bethe Ansatz states" with an 
arbitrary number of spins down. The relations of would then give the Bethe Ansatz 



equations. Some exact ground states can also be constructed using the theorem in [[L3f| . 



In the following we use a simpler method to solve the model. Our Hamiltonian (|6]) can 
be rewritten as (up to an irrelevant constant term and a constant factor): 

1*1 1 . 1*1 1 

H = - XX 2 + 2Sl '* ' ^l.i+lXg + 2S 2,« " S 2,i+l) - - XX 2 + 2Sl '* ' S 2,i+i)(2 + 2S 2.* ' S M+i) 
i=\ i=l 



5*1 1 

+ k £(« + 2S M • S 2 ,)(- + 2S lii+1 ■ S w ). 
6 J=1 2 2 



We define the rung states as 



(7) 



|0j >— — 7^(1 Tl,*j 4-2,i> — I Il,it2,i>), |lj >— | Tl,i) t2,i>? 



and the Hubbard operators = | ck, >< a, P — 0,1,2,3. Hamiltonian (0) can be 



|2i >— -7^(1 tl,i; j 2 ,i> +1 ll,i)T2,i>) ; |3j >— | | 2 ,i>) 



rewritten as 

1 N r on i i 



2 ^ 

where 



prr pdd , v00\/_ "yOO ^ 

r i,i+l r i,i+l + ^ A i ^i+l/ 



•^M+i — (2 + 2 ^i,« ' Si,t+i)(- + 2S 2 ,j • S 2 ,j+i), — (- + 2Si 5i • S 2 ,i+i)(- + 2S 2 ,j • Si ji+1 ). 

It can easily be checked that 

P% +1 \a t > >= \Pi > \a i+1 >, P% +1 \oh > lA+i >= (-l) £(a)+e(/3) |A > |a i+ i >, 

where e(a>) is the parity of the state \a >, e(0) = and e(l) = e(2) = e(3) = 1. The P 
operators can be expressed as 

prr +i = £ X?X&, PtUl = X (-lY (a)+m XfX^. 
a,/3=0 a,/3=0 
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Therefore, we can rewrite the Hamiltonian (6) or (7) as 

N 3 N -I 1 AT 

# = EEM + + - - - J £*?°> (9) 

i=l o=l i=l i=l 

with U = y> an d <^ = 0- 111 fact, the model (9) is integrable for arbitrary real constants {7 
and J. This corresponds to the case of an arbitrary rung exchange coupling and an arbitrary 
rung-rung biquadratic coupling. We shall discuss the generalized integrable case rather than 
the special case (7) in the following text. 

Obviously, Xf a represents the number of the local state |a, > on the rung i and satisfies 
the hard-core condition Y? a =o^t a = 1- ^ we choose |0 >= |0i > <8> • • • <8> |0tv > as the 
vacuum state, X"° (a = 1, 2, 3) can be looked upon the creation operators of a-particles, 
i.e., Xf°|0j >= |otj >. The particle numbers N a = Y^iLi Af a are conserved quantities. Notice 
that only three of them are independent since J2a=o N a — N. In this sense we construct the 
Bethe states \N±, N 2 , N3 >. As shown in Eq.(9), there is no hybridization among the states 
|1 >, 1 2 >, and |3 > because of the absence of Xp, Xf 3 , X/ 3 and their conjugates in Eq.(9). 
That means that the excitations from \Ni, N 2 , N 3 > to \N[, N! 3 > are dispersionless (i.e. 
have zero excitation energy) as long as Ni + N 2 + N 3 = N[ + N 2 + N 3 . In fact, the quantities 

N N 



-31 

=1 i=l 



commute with the Hamiltonian, which means that the Bethe states are highly degenerate. 
The general eigenstates can be constructed from \N e , 0, >: 

\N 1} N 2 , N 3 >= Y 2 fY.^\N ei 0, >, N e = N 1 + N 2 + N 3 . 
Therefore, we need only to consider the Bethe state \N e , 0, >= |7V e >. This state reads: 

\N e >= £ ^Kn^-.-^^jn^lO^ 

where ^(ni,n 2 , • • • , njvj is the wave function, and rij = 1, • • • , N denotes the coordinate of 
the j-th triplet rung. Let 77 be defined by U = 2 cosh 77. From an analysis similar to the one 



used in solving the XX Z spin chain (see, eg.|T4|, [15||), we have the Bethe Ansatz equation 
(for \ j ,j = l,..,N e ) 



N 



sin(A i — H sin(Aj - A; - vq) 

sin(Aj + §77) \ W sin(A i - Xi + irj) ' " 



and the eigenenergy to (0) (up to an irrelevant additive constant) 

2 sinh 2 r] 




(11) 



cosh rj — cos 2Xj 



where Xj are the rapidities of the triplet rungs. We note that a similar situation (namely 
mapping of a biquadratic spin-1 chain to an XXZ Heisenberg chain) was discussed in [16]. 

The phase diagram of the ground state spanned by J and U is almost the same to that 
of the XXZ Heisenberg spin chain with an effective magnetic field J, in the sense that the 
triplet rungs and the singlet rungs serve as the up spins and down spins, respectively. We 
distinguish three regions, according to U > 2, — 2 < U < 2, U < — 2 respectively: 

(i) U > 2: For \J\ < J c , the ground state is a Mott-like "insulator" consisting of N/2 
triplet rungs and N/2 singlet rungs with an energy gap (a gap at the lower end of the energy 
spectrum) A = J c — \J\, where J c is given by 



For J > U + 2, the triplet rungs are unfavorable and the ground state is a rung-dimerized 
state (product of N singlet rungs) with an energy gap A = J—(U +2), while for J < —(U +2), 
the ground state is a product of N triplet rungs with an energy gap A = | J\ — (U + 2). The 
latter two phases correspond to the completely polarized states in the XXZ spin chain. In 
the intermediate parameter region J c < | J| < U + 2 one has a gapless phase. 

(ii) —2 < U < 2: There is no Mott-like phase in this case. For J > U + 2, the ground state 
is still a rung-dimerized state (consisting of only singlet rungs) and for J < —(U + 2), the 
ground state is a product of iV triplet rungs. For |J| < U + 2, the ground state is a spin 
liquid with gapless spinon excitations (cf. [JTJJ for a discussion of phenomena of this type). 

(iii) U < —2: There is no gapless phase except for J = 0, U = —2. The ground state is 
almost the same as that of a ferromagnetic spin chain (in the same sense as above). For 
J < 0, the ground state is a triplet-rung product while for J > the ground state is a 
rung-dimerized state. 

Acknowledgements: We would like to thank Dr. R.H. Yue (North- West University, Xian) for 
very helpful discussions. The DFG and SFB-237 support to the second author is gratefully 
acknowledged. 



7r sinh rj 



V 



sech—(l + 2n). 

=-co ^T) 



(12) 
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